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Abstract 

We investigate multi-dimensional universe with nonlinear scalar curvature terms to 
evaluate the probability of creation of primordial black holes. For this we obtain Eu- 
clidean instanton solution in two different topologies: (a) S^^^ - topology which does 
not accommodate primordial black holes and (b) x S'^^^-topology which accom- 
modates a pair of black holes. The probability for quantum creation of an inflationary 
universe with a pair of black holes has been evaluated assuming a gravitational ac- 
tion which is described by a polynomial function of scalar curvature with or without 
a cosmological constant (A) using the framework of semiclassical approximation of 
Hartle-Hawking boundary conditions. We discuss here a class of new gravitational in- 
stantons solution in the i?*-theory which are relevant for cosmological model building. 
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1. INTRODUCTION : 

The last couple of decades witnessed a considerable progress in modern cosmology both 
theoretically [1] and observationally [2]. It has been generally accepted that inflation is one 
of the main criteria to build a consistent cosmological models of the early universe. The 
important outcome of inflation is that it explains satisfactorily the observed large scale 
structure of the universe in addition to solving longstanding problems both in particle 
physics and cosmology. The quantum fluctuation of the fields in the early universe may 
produce abundantly primordial black holes (in short PBH). In the literature the problem 
of creation and subsequent evolution of PBH in the early universe are discussed in different 
contexts [3]. In 1975, the discovery of Hawking radiation ushered in a new era in black hole 
physics [4] . The life times of black holes formed by gravitational collapse of a massive star 
are quite long which are in fact comparable to the age of the universe. However, primordial 
black holes are tiny and are comparatively short lived. Therefore, the hope of confirming 
Hawking radiation by observation is through PBH hunting. The existance of a pair of 
PBH may be important in the context of the recent cosmological observations. Bousso 
and Hawking (BH) [5] calculated the probability of the quantum creation of a universe 
with a pair of primordial black holes in a 4 dimensional universe considering two different 
Euclidean space-times : (i) a universe with space-like sections with S^- topology and (ii) a 
universe with space-like section with x S^- topology as is obtained in the Schwarschild- 
de Sitter solution. The first one describes an inflationary universe without a pair of PBH, 
while the later describes a Nariai universe [6] , an inflationary universe with a pair of black 
holes. BH considered a massive scalar field to investigate the problem which provided an 
effective cosmological constant for a while through a slow rolling potential (mass term) in 
the early epoch. However, the evolution of a universe at the Planck time t ~ Mp^ may be 
better understood in the framework of quantum gravity. A consistent quantum theory of 
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gravity is yet to emerge. It was Kaluza and Klein [7] who first independently initiated the 
study to formulate a unified theory of gravity and electro-magnetic interaction by intro- 
ducing an extra dimension. But the KK-approach did not work well subsequently because 
of uninteresting dimensional reduction technique employed in the theories. In recent years 
there is a paradigm shift in cosmological model building in the higher dimensions from that 
of previous approaches. The success of Superstring theory/M-theories revived the studies 
in higher dimensional framework which has been considerably generalized [8]. In modern 
higher dimensional scenario, the fields of the standard model are considered to be confined 
to (3-1-1) dimensional hyper-surface (referred to as 3-Brane) embedded in a higher dimen- 
sional space-time but the gravitational field may propagate through the bulk dimensions 
perpendicular to the brane which is referred to as braneworld [9]. At high enough energy 
scale, Einstein field equation in 4 dimensions may not give a correct predictions of the 
universe. Such an energy scale perhaps was available shortly after the big bang. Thus a 
more general theory is needed which at low energy reproduces the Einstein's theory. In the 
usual 4-dimensions or multidimensional universe the gravitational sector of the action may 
be modified by considering non- linear terms of scalar curvature {R). The motivation of 
studying higher dimensions is that most of the theories of particle interactions, including 
String theory requires space-time dimensions more than the usual four for their consistent 
formulation. It is therefore considered essential to check if consistent cosmological or astro- 
physical solutions, which can accommodate these theories are also allowed. In this paper 
we consider such a modified theory in multidimensional universe to evaluate the probability 
of pair creation of PBH using the prescription of Hartle-Hawking no boundary proposal 
[10]. According to the no boundary proposal, the quantum state of the universe is defined 
by path integrals over Euclidean metrics g^i, on compact manifolds M. From this it follows 
that the probability of finding a three metric hij on the spacelike surface dM is given by a 
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path integral over all g/^i, on M that agrees with hij on dM. With the inclusion of matter 
fields, one obtains the wave function of the universe which is given by 

where {hij, ^om) are the three-metric and matter fields on a spacelike boundary dM and 
the path integral is taken over all compact Euclidean metrics g^^. We now estimate ^ from 
a saddle point approximation to the path integral. Two types of topologies that describe 
inflationary universe characterized by spacelike sections, namely (i) R x S"^ -topology and 
(n) Rx X S'^ -topolo gy are considered. The former represents a universe without 
a pair of PBH and the later represents that with a pair of PBH. We obtain Euclidean 
solutions of the modified Einstein equation corresponding to each of these universes, which 
can be analytically continued to match a boundary dM of the appropriate topology. The 
Euclidean action / is evaluated for the above solutions. The wave function of the universe 
in the semiclassical approximation is given by 

^9m] ~ ^ne"^", (2) 
n 

where the sum is over the saddle points of the path integral, and denotes the corre- 
sponding Euclidean action. We can thus assign a probability measure to each type of 
universe: 

P[hij,^9M]^ei-''''''), (3) 

where the superscript Re denotes the real part of the action corresponding to the domi- 
nant saddle point, i.e., the classical solution satisfying the Hartle-Hawking (HH) boundary 
conditions [10]. 

It is well known that a theory with higher order Lagrangian is conformally equivalent to 
Einstein gravity with a matter sector containing a minimally coupled, self interacting scalar 
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field. However, the renormalisation of higher loop contributions introduces terms into the 
effective action that are higher than quadratic in R. Consequently it is important to study 
the effects of these terms in the quantum creation of a universe with a pair of PBH. Paul et 
al. [11] studied the probability of pair creation including term in four dimensions. Later 
an effective action upto term was considered in Ref. [12] and the probability of creation 
of a pair of black holes in the modified Einstein action in 4-dimensions is estimated. It 
is noted that the probabilty of creation of a universe with a pair of PBH suppresses in 
the i?^-theory without a cosmological constant. Subsequently, a higher dimensional uni- 
verse was considered in Einstein-Hilbert action to probe the PBH pair creation rate [13]. 
Giinther et al. [14] studied multidimensional gravitational models with scalar curvature 
non-linearities of types R~^ and R'^ to study the existance of at least one minimum of 
the effective potential for the volume moduli of the internal spaces with an emphasis to 
obtain stabilization. The stabilization of the internal space is important because the extra 
dimesional space components should be static or nearly static at least from the time of 
primordial nucleosynthesis. In contrast to i?^-models, the i?^-model shows a rich substrate 
of the stability region in parameter space which minimally depends on the total dimension 
of the bulk space. In this paper we consider a modified gravitational action described by 
a Lagrangian polynomial in R considering terms up to i?^ to probe the early universe. 

The paper is organised as follows : in section 2 we consider the action and the corre- 
spoding field equation is determined. We then obtain gravitational instanton solutions in 
two different topologies. In section 3 we use the action to estimate the relative probability 
of the two types of universes with different topologies in the presence of PBH. In section 4 
we give a brief discussion. 
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2. Gravitational Instantons with or without PBH : 



We consider a higher dimensional Euchdean gravitational action which is given by 



where g is the D-dimensional Euclidean metric, R is the Ricciscalar, A is the cosmological 
constant, and K = h"'^Kab is the trace of the second fundamental form of the boundary 
dM in the metric. We consider a polynomial function in Ricci scalar f{R) = Y^i^iR^, in 
particular Ao = -2A, Ai = 1, A2 = a, A3 = /?, A4 = 7, i.e., f{R) = R+aR'^+PR^+'yR'^-2A. 

A. S^~^- Topology, the de Sitter spacetimes: 
In this section, we study vacuum solutions of the field equation corresponding to the non- 
linear higher dimensional action (4). We look for a solution with space like section S^~^. 
We choose the D-dimensional metric ansatz, which is given by 



where d = D — 1, ais the scale factor of a D dimensional universe and dCl'^ is a line element 
of unit {D — l)-sphere. The scalar curvature is given by 



where the overdots denotes differentiation with respect to r. We use the constraint through 
a Lagrangian multiplier (3 and rewrite the action as 





dS'^ = dT'^ + a^{T)dnl, 



(5) 




(6) 




/ d^^'-^^xVhf'iR) 

JSM 



(7) 



where Vg 



1 27r(''+l)/^ 
167r r(d+l)/2 



. Variaton of the action with respect to R, determines (3, which is 



P = a''f'{R). 



(8) 
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Substituting the above constraint in the action we get 



Ie 



= -Vo / 

Jt=0 



d2 1 



f{R) - f'{R) [R-d{d-l)[^ + ^]\+ 2dRf"{R) 



dr + 



2dVoaa-'-'f{R) . (9) 

J T=0 



Varying the above action with a and R respectively, we get 

1 



f'iR) 
+ f"{R) 



2(d-l)- + (d-l)(d-2) , „ „ 



{d-l)-R + 2R 



f /// r T->\ 62 



+ R 

+ 2f"{R)R' - f{R) = 



a 

2d-+d{d-l) (^-^]+R 



0. 



(10) 
(11) 



The above field equations admits an instanton solution given by 



a = — sin (i?oT), 



where Hq is determined from the constraint equation 



(12) 



d{d-l)Hl+ad^{d+l){d-2,)H^+(3d\d+lf{d-h)Hl+-id*{d+lf{d-7)Hl-2K = 0. (13) 

Thus Ho = f{d,a,P,j,A) determined by five parameters. We note that for (i) D = 
4, Ho = /(c?,/?,7,A); (ii) D = 6, Ho = /(d,a,7,A); and (iii) D = 8, Ho = f{d,a,l3,A). 
The above instanton solutions satisfy the HH boundary condition viz., a(0) = 0, d(0) = 0. 
Now one can choose a path along the r^*' axis to r = the solutions describe half of the 
Euclidean de Sitter instanton in -topology. Analytic continuation of the metic (5) to 

Lorentzian region xi ^ + ia gives 



ds'^ = dr'^ + a^(r) — d<T^ + cosh^a dQ^_2 



(14) 



which is a de Sitter like metric. However, if one sets t = it and cr = ^ + X) the metric 
becomes 

ds'^ = -df + S^{t) \dx^ + sinh^x d^l-2 \ > (15) 
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Re _ 



where S{t) = —ia{it). The above metric now describes an open universe. Thus, the creation 
of an open inflationary universe may be reaUzed in this case. It may be pointed out here that 
a closed universe created quantum mechanically may be realized as an open inflationary 
universe under a Wick rotation. Since it is not yet known whether our universe is exactly 
flat, it may be interesting to study an open universe. The real part of the Euclidean action 
corresponding to the solution calculated by following the complex contour of r suggested 
by BH is considered here which determines 

'd{d + l)Hl + a(f{d + ifH^ + Pd^{d + if Hi + 7d^(d + if Hi - 2A] , 

(16) 

where Id = Jq"° sin'^j/ dy, denoting y = HqT. The value of the integral Id depends on the 
dimensions of the universe. For odd number of dimensions (D), 

^ d-ld-Z 3l7r ^ ^ ,^ . 

id= — ; ; — - TTTTT, for a even, (17) 

d d-2 422' ' ^ ^ 

and for even number of dimensions (D), 

^^=^^ M<io<i.. (18) 

With the choosen path for r, the solution describes half the de Sitter instanton in a 
higher dimensional universe with S'^ topology, joined to a real Lorentzian hyperboloid 
of {R^ X S^) -topology. It can be joined to any boundary satisfying the condition agM > 0. 
For aQM > Hq^^ the wave function oscillates and predicts a classical space-time. 
B. 5^X5^-2- Topology : 

We consider vacuum solution of the field equation corresponding to the action (4) and 
look for a universe with S^XS^~^ -spacelike sections. The above topology accomodates a 
pair of black holes. The corresponding metric ansatz for D = 1 + 1 + d dimensions is given 

by 

ds'^ = dT^ + a^{T)dx^ + b^{T)dnl, (19) 
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where a(r) is the scale factor of two sphere and 6(r) is the scalefactor of the ^ ^-sphere 
given by the metric 



= dxi + siv?x\dx\ + siv?' xisiv?' X2dx\ + .... d space. 



The scalar curvature is given by 



R 



^ 'd ^ ,b , f b'^ 1 \ ^ ,db 



(20) 



Using the above constraint (20) the Eucledian action (4) may be rewritten as 



Ie 



-V 



f{R)ab'' -p(R + 2U2dl + d{d 2d^^ 



dr 



OTT JSM 



(21) 



where = r((D+iy2) ■ undetermined multiplier ^ may now be obtained by varying 
the action and inserting the corresponding (3 in the action, we finally obtain 



fe2 



,db 



1 



Ie = -V;' / ab'^i f{R) - f'{R) [R-d{d-l)r^-2d-- d{d - 1)- ) ) dr - 



r=0 



'62 



d ,6 
a b 



■ (22) 



r=0 



Varying the action with respect to a, b and R respectively, we get the following field 
equations : 



f"{R) 



R + 2^ + 2dl + 2d^+d{d-l) ( ^2 

a b ab V o o , 



0, 



(23) 



f'{R) 



b (i? 



+ 2f"{R) [R + 2dR'^ + 
2R?f"\R) - f{R) = 0, 



(24) 



f'iR) 



i. + 2(d-l)^ + 2(d-l)^ + (d-l)(d-2)a-^U2^ 



+ 



2f"'{R)R^ + 2f"{R) 



(d-l)Rj+R + R- 
b a 



- f{R) = 0. (25) 
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We now look for instanton solutions from the above field equations which is given by 

a = 4 sm{HT), b = H'^, R={2 + d)H^ (26) 

H 

where H satisfies the following constraint equation 

dH'^ + a{2 + d){d - 2)H^ + /3(2 + df{d - 4)H^ + 7(2 + df{d - 6)H^ - 2A = 0. (27) 

Here H depends on the parameters a,/?, 7, A i.e., H = /(d, a, /?, 7, A) for a given extra 
space dimension (d). We note that for (i) D = 4, H = /(c?, /?, 7, A); (ii) D = 6, H = 
/(d, a,7,A); and (iii) D = 8, H = f{d,a,P,A). The above solution satisfies the HH 
boundary conditions a(0) = 0, d(0) = 0, 6(0) = 0, 6(0) = 0. Analytic continuation of 
metric (19) to Lorentzian region i.e., using Wick rotations t ^ it and a; — | + icr, one gets 

ds'^ = -dt^ + S{tf da^ + H'^ dnl_2, (28) 

where S{t) = —ia{it). In this case it is evident that the analytic continuation of time and 
space is not an open universe. The space-time now represents an anisotropic universe. The 
corresponding Lorentzian solution is given by 

a(Ti™) Lrb-^ = H-^coshHT^"^ 
6(ri"^)|^Re-^ = H-^ 

Its spacelike section represents {D — 2)-spheres of radius a with a hole of radius 6(= H~^) 
punched through the north and south poles. The physical interpretation of the solution 
is that of {d — 1) spheres containing two blackholes, they accelerate away from each other 
with the expansion of the universe, at the opposite ends. The real part of the action may 
now be determined following the contour approach suggested by BH, which is 



((2 + d)H^ + «(2 + dfE"" + f3{2 + dfH^ + 7(2 + d)^H'' - 2A) 

(29) 
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where d = D — 2. The solution (26) describes a universe with two black holes at the poles 
of a (Z) — 2)-sphere. It may be pointed out here that the non-zero contribution to the 
action comes from the surface terms only. 

3. Evaluation of the probability for PBH : 



In the previous section we compute the action for inflationary universe with or without a 
pair of black holes. The probability for creation of a higher dimensional de Sitter universe in 
/(i?)-theory may now be obtained using gravitational instantons obtained in the previous 
sections. The probability for nucleation of a higher dimensional universe without PBH is 



\D[D-l)Hl+aD'\D-lf H^+l3D'-i{D-lf H'^+'yD^{D-l)^Hl-2A\ 



(30) 



where Ho satifies the constraint eq.(13). However, the probability of nucleation of an 
inflationary universe with a pair of black holes is obtained from action (22), which is 



51x50-2 



2V<;(D-3)(0-2)/2 

TT^ 



(31) 



where H satifies constraint eq.(27). 



Special Cases : For D > 4 we note the following : 
(i) a = /3 = 7 = 0, the probability estimates given in (30) and (31) reduce to 



PgD-i ~ e 



4^,^,^(£-^iD-2)^(0-2)A 



(32) 
(33) 



We recover the probabilities obtained by Bousso and Hawking in Einstein theory when 
D = A, which is 



Pgs ~ e^''/^, P51X52 ~ e^""/^. (34) 
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In the above de Sitter universe is more probable for a positive cosmological constant. 

(ii) For a = (3 = 0, 7 / 0, the corresponding probabilities are 



SB- 



-1 ~ e 



Vq'd^i ( -&d(d-i)hI 

0-8 



+16A ^ 



51x50-2 ~ e 



IT^ / 



(35) 
(36) 



where Ho and H are determined from eqs. (13) and (27) respectively for D 8. In four 
dimensions the probabilities reduce to 



Pg3 ^ e 



(37) 



For A = and D 8, we evaluate the probabilities which are 



PcD-i ~ e^^^yolD-i{D{D-l))'^/^{8-D)('^-^)/^{^)('^-^)/<^] 



D ri2V'DO/2(8-ri)(0-8)/6m-3)('f'-2)/2(-X-){0-2)/61 



(38) 



(39) 



We note the following: 

(i) a new gravitational instanton solution in higher dimensions is obtained for a i?^-theory 
with a negative couphng parameter ( i.e., 7 < 0) for D > 8, and that with a positive (i.e., 
7 > ) for D < 8, (ii) de Sitter universe is more probable for D > 3. 

We also note that for = 8, an instanton solution with = ^ in S^-topology and with 
H"^ = in 5^ X S'^-topology are permitted. The probabilties of PBH for D = 8 dimensions 



are 



Pg7 ~ e 



if||5-(27+20487A3) 



(27+20487A3)j 



(40) 



We further note that the probabilities are same if (i) I7I = 2048A3 A > with negative 
7 and (ii) 7 = 2048A^ A < with positive 7. An interesting gravitational instanton is 
permitted without a cosmological constant in D = 8. 
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(iii) For a 7^ 0, /3 = 7 = 0, the probabilities are 



50- 



-1 ~ e 



51 



-2 ^ e 



!Vo/i3_i /-2D{D-l)ff2+8A\ 

7 

2V„'(i3-3)(g-2)/V-20if2+8A \ 

) 



(41) 

(42) 



In this case PBH pair creation may be possible even without a cosmological constant. For 
A = 0, the probabilities are given by 



so- 



-1 ~ e 



(D-4)/2 



4yj(§5|)'°-'''/''i30/2(4-i3)(-D-4)/2 



(43) 
(44) 



We note that the instanton solutions considered here for the two different topologies phys- 
ically sensible when a < for D > A and a > for D < A. We note that in higher 
dimensions with space-time dimensions D > 3 the probability of nucleation of a S^~^- 
topology is found to be more than that of a universe with x S^'"^ topology. 
In 4-dimensions the probabilities reduce to 



P53 ~ e 



[247ra+37r/A] 



[167ra+27r/A] 



(45) 



A theory with a positive cosmological constant and a > 0, leads to a universe without 
PBH, which is more probable. Though a negative ct < — ^ leads to greater probability for 
a universe with PBH, the negative values of a lead to a classical instability in i?^-theory 

[11]. 

(iv) For a = 7 = 0, /3 7^ 0, the probabilities are 



50-1 ~ e 



Six50-2 ~ e 



Volp-i f -iD{D- 



l) wg+12A \ 



2V^(D-3)(-P~2)/2 



(46) 
(47) 
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with D 7^ 6. The probabihties in D = 4 dimensions may now be obtained as 



(48) 



We note that PBH pair creation may be described even without a cosmological constant. 
The corresponding probabilities are 



P50-1 ~ e 



8K,/D_i(-^)<''-'^/^D(i?-l))0/2(6-i3){0-6)/4] 



8V; 



;'(i:>-3)(-D-2)/2(-g^)'°-''/*i30/2(6_i:,){i3-6)/4 



(49) 



(50) 



We note that a universe with space-time dimensions (i) D > 6 admits an instanton solution 
with 7 < and (ii) D < 6, with a positive j3. The probability measure for pair creation in 
inflationary universe is more comparable to that with a universe without primordial black 
holes. 

Let us consider D = 6 dimensional universe. In this case = ^ in 5^-topology and 
= ^ in 5^ X S'^-topology. The corresponding probabilities are 



Pg5 ~ e 



■(1+27/3A2 



^(l+27^A2) 



(51) 



de Sitter universe without a pair of PBH is more probable. 

(v) For A = 0, the probability measure in the two topologies under consideration are given 



by 



P5C-1 ~ e 



2VoJp-l 



2aD'^(D-l)^H^+4pD^(D-l)^H2+6jD'''(D-l)^H^ 



2Vo'(D-3)(-''-2)/2 

ir^ 



2oiD'^H*+4,Pd'-^H^+G-,D^H^ 
D-2 



(52) 
(53) 



In D = 4 dimensions the probabilities reduce to 



P53 ~ e[''(^^'^+^^^^-^°+^°^^^^-^°^],P5ix52 ~ f,Hie'^+i2»fiH^+76S^H^)]_ 



(54) 
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We now discuss few other special cases which are interesting and follow from the general 
result obtained above in higher dimensions. 



(a) For 7 = 0, the probabilities are 



P50-1 ~ e 



2c,d'^(D-1)'^H^+4I3d'-^(D-1)^ 



2V'(i3-3)(^-2)/2 

IT^ 



D-2 



(55) 
(56) 



one determines the instanton solution even with A = it reduce to that obtained in Ref. 
[9] for L> = 4 

(b) For P = 0, and A = 0, the probabilities are 



— TTD — 



2aD'^(D-lfH^+6-yD*(D-lfH^ 



P. 



2Vj(D-3)(0-2)/2 

IT^ 



2aD'^H*+6-tD^H^ 



For £> = 4 



P53 ~ e[3'^(8"+3(^287)i/3)]^^^^^^^ ^ g[27r(8a+3(1287)i/3)]_ 



(57) 

(58) 
(59) 



Here = ^/^^S^ ~ '\/t^^ solutions. In this case a physically interesting 

instanton solution is obtained for a positive value of 7. For a < and if |8a| < 8(1287)^/^, 
de Sitter universe without PBH is more probable. However interesting possibilities emerges 
when |8a| > 8(1287) V3. 



4. Discussions : 

In this work we studied the probability of pair creation for primordial black holes in a non 
linear theory of gravity. We consider an action with a polynomial function of Ricci scalar 
R as f{R) = YLi^iR^ to obtain gravitational instanton solutions in the multidimensional 
universes. The Euclidean action is then evaluated using instanton solutions in the two 
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cases: (i) a universe with R x S'^-topology and (ii) a universe with R x x ^'''-topology 
in higher dimensions. The former admits an inflationary universe without PBH and later 
a pair of PBH. We present the corresponding estimation of the probability in a multidi- 
mensional universe and compare then to predict the probability of creation of universe. 
We determine that the probability measures in both the topologies in the framework of 
a non-linear theory of gravity with square, cubic and quadratic power of scalar curvature 
term (R) in the gravitational actions and estimate the coupling parameters for a physically 
realistic instanton solutions. One interesting point is that in this f{R) gravity model with 
square and higher order curvature terms one obtains gravitational instanton solution even 
without a cosmological constant. We note that the de sitter universe without a pair of 
black holes is more probable for D > 3. We found that the probability of a universe with 
topology R X S^- turns out to be much lower than a universe with topology Rx x in 
ii!^-theory unless A < — ^ in 4-dimensions. The probability of creation of a universe with 
a pair of PBH is found to be suppressed if the R^ term is included in the action under the 
constraints |a| < 2y/]3 or a > —2y/]3 with A = 0. Again the probability of creation of a 
universe with a pair of PBH is found to be strongly suppressed if one extends the polyno- 
mial function of f{R) gravity model upto i?^-term with the constraints |8a| < 8(1287)^/^ 
or 8a > —3(1287)^/^. One obtains Hawking- Turok [15,16] type open inflationary universe 
in R X S^-topology, but in the other topology accommodating PBH does not permit an 
open universe instead it gives rise to an anisotropic universe. We note a dimensional de- 
pendence of the gravitational instanton and deterimned the different parameters to realize 
such instantons in the non-linear theory of gravity. A class of new gravitational instantons 
are presented here which are relevant for cosmological model building. 
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